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ABSTRACT 

We present estimates of the energy input from supernovae (SNe) into the intergalactic 
medium using (i) recent measurements of Si and Fe abundances in the intracluster 
medium (ICM) and (ii) self-consistent gasdynamical simulations that include processes 
of cooling, star formation, SNe feedback, and a multi-phase model of the interstellar 
medium. We estimate the energy input from observed abundances using two different 
assumptions: (i) spatial uniformity of metal abundances in the ICM and (ii) radial 
abundance gradients. We show that these two cases lead to energy input estimates 
which are different by an order of magnitude, highlighting a need for observational 
data on large-scale abundance gradients in clusters. Our analysis indicates that the 
SNe energy input can be important for heating of the entire ICM (providing energy 
of 1 keV per particle) only if the ICM abundances are uniform and the efficiency 
of gas heating by SN explosions is close to 100% [esN ~ 1, implying that all of the 
initial kinetic energy of the explosion goes into heating of the ICM) . 

The SNe energy input estimate made using simulations of galaxy formation is con- 
sistent with the above results derived from observed abundances, provided large-scale 
radial abundance gradients exist in clusters. For the cluster AWM7, in which such a 
gradient has been observed, the energy input estimated using observed metal abun- 
dances is ~ 0.01 and ~ 0.1 keV per particle for esN — 0.1 and esN = 1, respectively. 
These estimates fall far short of the required energy injection of ~ 0.5 — 3 keV per 
particle that appears to be needed to bring models of cluster formation into agreement 
with observations. Therefore, our results indicate that, unless the most favorable con- 
ditions are met, SNe alone are unlikely to provide sufficient energy input and need to 
be supplemented or even substituted by some other heating process(es). 
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1 INTRODUCTION 

Hierarchical models of structure formation have been very 
successful in explaining many observed properties of galaxies 



and galaxy clusLeis. Never Llieless, some puzzling piublems 
DLULial LliLuiLLital jLudicj liauL dLmuniLiiilLd 



remam upcn 
that so me heating of gao 



in addition to the heating dur 



ing the gravitational collapse, is required to explain t he ob- 
served properties of the intracluster medium (ICM). Kaiser 



(1991) pisL showed LliaL an eaily injeLLlun uf energy lesuUs 



in corr olationa and evolution of bulk oluatcr propcrtica (X 
ray luminosity, gas temperature, etc.) that match observa- 
tions. This con clusion was backed by numerical simulations: 
Evrard (1990) was able to get a better fit to the X-ray lu- 
minosity of the Coma cluster in his cosmological si niulation 
by preheating the gas to 10^ K (~ 0.9 keV), while Navarro 



□h 



ubble Fellow 



et al. (1995) showed that pre-heated clusters matched the 
observed slope of the correlation between X-ray luminosity 
and temperature. Recent semi-analytical st udies of cluster 
evolution have reached si milar conclusions (|Cavalicrc et al 



1997; Balogh et al. 1999a; Valageas & Silk 1999 



Wu et al 



T999| ) 

The exact amount of required energy injection depends 
on the epoch, and have be en argued to be in the range of 



Cavaliere et 



0.5 — 3 k e V per gas particle ([Navarro et al. 1995 
al. 1997|; |Balogh et al. 1999a|; [Wu et al. 1999^ Ahhough the 



problem has been identified, it is not yet clear what processes 
can provide the required heating. It is clear that identifica- 
tion of these processes is crucial for a complete picture of 
cluster formation. 

The candidate process which has been discussed most 
is supernovae-driven galactic winds. The gas of the galac- 
tic interstellar medium (ISM) can be heated by supernovae 
explosions and acquire energy comparable to or larger than 
its gravitational binding energy. This heated gas can then 
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flow away and result in additional heating if the winds result 
in shocks when they encounter intergalactic gas. Although 
there is observati onal evidence for suc h winds in present- 
day galaxies (e.g., Heckman et al. 199C), theoretical models 
of winds are rather ill-constrained due to uncertainties in 
the efficiency of conversion of supernovae (SNe) explosion 
energy into thermal energy of the gas and other details. 

Early estimates of possible energy input from SNe based 
on observed metal abundances showed that SNe are plausi- 



of Anders fc Grevesse (1989) {npe/uH = 4.68 x lO"'^ and 



ble can didates (e.g., White 1991 



David et al. 1991 



Loewen- 



stein & Mushotzky 1996). However, in these estimates it 



was assumed that the distribution of metals in the ICM is 
uniform and that the efficiency with which energy of SNe 
explosions can be converted into thermal energy of the gas 
is close to 100%. Therefore, there is a need for detailed es- 
timates using new measurements of metal abundances in 
clusters and current galaxy formation models which have 
become much more advanced and sophisticated in the last 
several years. 

Evaluations of SNe as a heating source hav e recently 



been performed by Valageas fc Silk (1999) and Wu et al 



(1999)[]ising a semi-analytical approach to galaxy modelling 
These authors conclude that i t is unlikely that super novae 
are the only source of heating. Valageas & Silk (1999) argue 



then that radiation from quasar population can provide the 
required heating much more easily. 

In this paper we repeat previous estimates of the pos- 
sible SNe energy input using updated values of observed 



ICM at undances and relaxing the assumption of abundance 



unif orrmty, motivated by rec ent observations ( hizawa ct al 
1997; Finoguenov et al. 200C). We also make a separate esti- 
mate for the cluster AWM7, for which the radial abundance 
gradient was measured. The details of this estimate are pre- 
sented in § ^. We complement this analysis with direct es- 
timates of the energy input by counting the total number 



of supernovae exploded in all cluster galaxies throughout 
their evolution in self-consistent three-dimensional gasdy- 
namical simulations of galaxy formation. The simulations 
include cooling, star formation, SNe feedback and a multi- 
phase model of the interstellar medium in galaxies and have 
been shown to match many fundamental observed corre- 
lations of galactic properties such as the galaxy luminos- 
ity function, the TuUy-Fisher relation and its scatter, the 
color-magnitude sequence, and, perhaps most importantly, 
the evolution of the global star formation rate in the Uni- 
verse. The details of the simulations are described in § ^ 
The energy input estimates are presented and compared in 
§ ^ and discussed in § ^ We summarize our main results and 
conclusions in § |^ 



2 SUPERNOVAE ENERGY INPUT FROM 
OBSERVED ICM METALLICITIES 

We will first estimate the energy input from SNe to the 
ICM using observed metallicities of the cluster gas. We have 
based the estimate on silicon (Si) and iron (Fe) abundances 
because these two elements have been mos t accurately mea- 



sured for a large sample of g alaxy clusters ( Mushotzky et al 



1996; Fukazawa et al._199£ ). We use average ICM metallic- 
ities quoted in Table |l| and photospheric solar abundances 



nsi/uH = 3.55 X 10-^). 

Given that the mass of an element Xi, Mjjf. , within the 
cluster virial radius is known, the number of SNe type I and 
II required to produce this mass is equal to fx-Mxi/yi{Xi) 
and (1 — fx )Mxilyii{Xi), respectively. Here fx. is mass 
the fraction of the element contributed by type I SNe, yi {Xi) 
and yii{Xi) are the mass- weighted yields of the element Xi 
by SNe type I and II respectively. The SNe energy input can 
then be obtained by multiplying the number of SNe by the 
energy transferred to the gas during each SN explosion: 



Xi 



^SNlEsN 



,(1) 



where we denote Esni and Esnii energies released in ex- 
plosion of the two types of SNe, and esNi and esNii are 
fractions of the released energy left after the radiative losses 
during and after the explosion, which can be actually trans- 
ferred in the form of thermal and kinetic energy to the ambi- 
ent gas and lead subsequently to the increase of its entropy. 
There is a varying degree of uncertainty in our knowledge 
of the above parameters. 

First of all, our estimate of the mass of an element 
depends on the assumption about uniformity of the ob- 
served metallicities. If strong radial abundance gradients ex- 
ist in clusters, the observed metallicity is emission-weighted 
and therefore corresponds to the metallicity in the cluster 
core. Numerica l simulations of cluster formation (Metzler 



& Evrard 1994; Metzler & Evrard 1997) that include mod- 



elling of galaxy feedback predict the existence of strong ra- 
dial metallicity gradients in the ICM, as well as patchy spa- 
tial distribution of metals. At present, however, it is not clear 
whether large-scale abundance gradients are universal in 
clusters. Although abundan ce gradients have been observed 
in several cl usters (see, e.g., Allen fc Fabian 199^ ; Dupke fc 
and references therein), these are usually clus- 



White 1999 



ters that have a central cD galaxy and exhibit signatures of 
a central cooling flow (Allen & Fabian 1995). The spatial 



extent of the observed gradients coincides with that of the 
cooling flow region. It is thus unclear whether such central 
gradients imply the existence of a larger-scale gradient or 
they are simply due to the presence of a central cD galaxy 
and cooling flow. 

Currently, abundances in the fainter, outer parts of clus- 
ters can be me asured only for brigh t nearby systems. In 



a recent study, Ezawa et al. (1997), found strong large- 



scale metallicity gradients in the nearby cluster AWM7. 
The observed iron abundance in this cluster decreases from 
^ 0.5±0.05 solar within the central 60h~^ kpc to « 0.2±0.2 
at 300 — 500h~^ kpc. The radially averaged gradient can be 
well fitted by a /3-model with a core radius equal to that of 
the gas and /3 — 0.8. The Ezawa et al. measurement was the 
first in which the abundance gradient has been found far be- 
yond the cluster core radius. I t is not yet clear how com mon 
such large-scale gradients are. Pinoguenov et al. (2000) show 



that large-scale metallicity gradient are indeed observed in 
many clusters. It is clear, however, that strong large-scale 
gradients are not universal; for exam ple, no strong grad i- 
ent was detected in the Coma cluster ( Hughes et al. 1993 ). 
Clusters may therefore exhibit a variety of metal distribu- 
tions a nd span a range in the I CM metallicities. In support 
of this, Allen fc Fabian (1998) , present evidence that clus- 
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ters without strong metallicity gradients have systematically 
lower metal abundances. 

For our purposes it suffices to consider two possible ex- 
treme assumptions about the metal distribution. In real- 
ity the mass of metals will likely lie in between the masses 
computed under these assumptions. The first assumption is 
that the metallicity of the ICM is spatially uniform. Ob- 
servationally, the metallicity derived from a spatially unre- 
solved spectrum is emission- weighted. It is clear then that if 
a strong metallicity gradient is present in a cluster, the to- 
tal mass of metals may be significantly overestimated under 
the assumption of spatial uniformity. Our second assump- 
tion is that the metallicity gradient of the form observed 
by Ezawa et al.: Z{r) = Zo[l + {r /r^f]'^'^^^ is a universal 
property of th e cluster ICM. We w ill assume rc = 100/i~^kpc 
and P = 0.8 jEzawa et al. 1997| ) for all clusters, normaliz- 
ing Zo to a value such that Z{rc) is equal to the observed 
value of metallicity. Most of the cluster emission comes from 
radii < 2rc, providing an approximate way to account for 
the emission-weighting of the metallicity. The core radius 
of 100/i~^ kpc is larger than the best fit value for the clus- 
ter AWM7 but is closer to a typical core radius of the gas 
distribution in rich clusters. In addition to the estimate for 
the whole range of cluster masses, we will present the SN 
energy input for the specific case of AWM7 for which the 



abundance gradient has been observed and its parameters 
measur ;d. 



Another source of uncertainty is the relative importance 
of type la SNe (parameter fx, ) in the meta l enrichment of 



the ICM (Loewenstein fc Mushotzky 1996| ; Gibson et al 



1997; Nagataki & Sato 1998). Therefore, in the case of iron, 
we will treat fp^ as a free parameter and calculate Esn for 
values {fx = 0.0, 0.5, 1.0). Silicon is a special case, because 
yii{Si) ~ yi{Si). This renders Esn almost insensitive to a 
particular choice of This insensitivity, together with the 
fac t that silicon abundance was fairly accurately measu red 



by Mushotzky et al. (1996) and Fukazawa et al. (1998) 



ef- 
fectively reduces the uncertainties and thus makes Si a very 

useful plpmpnt fnr nnr patimat.P 



Th e third major aourcc of uncertainty ia the yiclda 



dieted by differe nt theoretical models of SN explosions ( Gib- 
son et ,J. 1997). The yields of SNe type II may depend on 



the initial metallicity of SNe, input physics of a model, and 
ot her factors. In our estimat e we will use yields calculated 
by [Woosley fc Weaver (1995)| for metal-poor {Z/Zq = lO""*) 
SNe with explosion energy of « 1.2 x lO^^ergs (model A) 
and metal-rich [Z/Zq = 1) SNe with explosion energy of 
« 1.2 X lO^^ergs for SN of mass < 25M0 and « 2 x lO^^ergs 
for SN of mass > 25M0 (Model B) . Model B has a higher ex- 
plosion energy for very massive stars to reduce the effects of 
reimplosion of explosively synthesized ejecta, thus increas- 
ing the yields. The yields for these models approximately 
give the lower and upper limits of the current theoretical 
predict ions (see Gibson et al. 1997 ; Nomoto et al. 1997a| ; Na- 



gataki ^ Sato 1998 ). Given the wide spread in predictions 



of theoretical models, the supernovae energy input estimates 
from the observed metallicities have the uncertainty of up 
to ~ 50%, in addition to other possible uncertainties. 

We calculate the average yield of SNII by averaging the 



mass-dependent yields with a stellar initial mass function 
(IMF): 



X^" yii{Xi,m.)(j}{m)dm 
II" 4>(m)dm 



(2) 



For the IMF we assume the (Salpeter 1955) function, 
4>{m) oc ra"^ "^^, with the lower mass limit of mi — 12Mq 
and mi — IIMq for the two yield model^ and the upper 
mass limit of mu ~ 40Mq. The SN with masses < IIMq do 
not contribute significantly to the metal enrichment, while 
stars of mass > 40Mq are rare. Our choice of the Salpeter 
IMF does not affect the average yields significantly: averag- 
ing with considerably shallower {(j>{m) oc m~^) and steeper 



(0(m) 



IMFs results in average yields that differ 



by less than 10% from the values for the Salpeter IMF. For 
SN type la, the yields appear to be independent of the SN 
mass. We use SNIa yields (see T able 1) predicted by the W7 
model of Nomoto et al. (1997b), which is a model of simple 
deflagration. The yields for the supernovae of type la and 
type II used in our analysis are summarized in Table ^. 

The energy released in a SN explosion {Esni and 
EsNii) also depends on a variety of factors (e.g. the mass 
of supernova). However, it can vary only by a factor of ~ 2 
and we will therefore assum e for simplicity that Es ni ~ 
EsNii = 1.2 X lO'^^ergs (e.g., [Woosley fc Weaver 1986|). Not 



all of this initial kinetic energy of expl osion is retained by 
the ejected g as. Analytical arguments ( Larson 1974 ; Babul 
& Rees 1992) and recent numerical sirnulations ( [Thornton 
et al. 1998p suggest that at most ~ 10% of the initial ki- 



netic energy of explosion can be ultimately transferre d to 
the ambient gas. In particular, Thornton et al. (1998) have 



numerically studied the evolution of the ejected material for 
a variety of densities and metallicities of the ambient gas and 
concluded that regardless of the ambient density and metal- 
licity, > 90% of the initial energy acquired by ejecta during 
the explosion is lost to radiation (see, however, discussion 
in § l^l). Therefore, we will assume that only 10% of the ex- 
plosion can actually be transferred to the ambient gas: i.e., 
esjv/ = ^SNii ~ 0.1. Here, we neglect the dependence of e on 
environment, metallicity, and possible systematic differences 
between esN i and esNii- However, according to Thornton 
et al. (1998), radiation losses are ~ 90% or more for most 
of the realistic environments and metallicities and the value 
of e = 0.1 is a reasonable upper limit for both types of SNe. 
We note that the energy estimates from observed metallici- 
ties presented in Figs, km can be simply linearly scaled up 
or down for other values of e. Clearly, the energy of SNe 
explosions is actually released into the interstellar medium 
and needs then to be somehow transferred to the IGM. Even 
if such transfer is possible, it is likely that it would result 
in additional energy losses. The estimates we make should 
therefore be considered as the upper limits on the amount 
of the energy that could have been available for the IGM 
heating. 

All of the estimates are made assuming the low-density 
flat cold dark matter model with cosmological constant 
(ACDM). The contributions of baryons, cold dark mat- 



t We chons p nnt to iisr; extrannlation and use mass limits of the 
yield grid of Woosley fc \Ycavcr f l9951 . This d oes not significantly 



affect the average y ields. Gibson et al. (1997)| and 



Loewenstein & 



Mushotzky (1996) have used a somewhat larger range of masses 



and obtained similar average yields. 
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Table 1. Parameters adopted in the estimate of energy input 
from observed metallicities 





Z 


Vl 






Si 


0.65 


0.158 


0.124 


0.158 


Fe 


0.32 


0.744 


0.096 


0.153 



reproduce many of the observed galactic properties which 
provides support to the used A'^sn — Mb relation. 

The simulations of the COB_E-normalized ACDM model 
(fio = 1 - Ha = 0.35; D.b = 0.026; Ho = 70 km s"^ Mpc'S 
where flo, fib, and Ho are present day values of the mat- 
ter and baryon densities and the Hubble constant, respe c- 



ter, and vacuum energy are: Qt ~ 0.05, flm ~ 0.25, 
Qa = 0.7, respectively. We assume a Hubble constant of 
Ho = 70 km s~^ Mpc~^. The cluster virial radius and mass 
for this model are defined at the overdensity of ~ 334 and 
we assume the baryon fraction within the virial radius is 
fb = ^b/Qm « 0.17'. 



tively) used here are described in Elizondo et al. (1999b) 
A total of 11 simulations were run from different realiza- 
tions of initial conditions. The size of the simulation boxes 
was fixed to Ltox = 3.5/i~^ Mpc = 5Mpc and the simu- 
lations were run using 128'^ grid cells and particles which 
gives mass and spatial resolution of ~ 2 x 10^ Mq and 
~ 27h~^ kpc, respectively. A total of 240 galaxies, 140 of 
which have Mb{z = 0) < —14, were formed in all the runs 
combined. The observed color-magnitude diagram, luminos- 
ity function (LF), and TuUy-Fisher relation of low-redshift 



galaxies are reproduced well by the simu lated galaxies (Eli 



3 F.MFRnV mPTIT FROM STTPFBIVnVAF, TINT 



NUMERICAL SIMULATIONS OF GALAXY 
FORMATION 

The question we now ask is what energy, Egj^, can be ex- 
pected to have been released in all the SNe in cluster galax- 
ies throughout their evolution? This question can be an- 
swered only in the framework of a self-consistent model of 
galaxy formation. There are currently two independently 
developing approaches to modelling galaxy formation and 



zondo et al. 1999c ; Elizondo et al. 1999b ). Simulations used 
here were done assurn ing SN feedback parameter of ^ = 200 



(see Yepes et al. 1997 ). This value of the parameter means a 
moderate efficiency of supernovae feedback and, correspond- 
ingly, relatively high star formation rate. 

The redshift dependence of the global star formation 
rate averaged over all simulations is shown in Fig. |^ together 
with the current data on obse rved global star fo rmation his- 
tory in the Universe (see also Yepes et al. 1999 for compar- 



evolution: semi-analytic mode l s (SAMs; e.g., [ Kauffmann et 
al. 1993| ; |Baugh et al. 1996[ |Somerville fc PrimackH^ 
and r eij orcnoca therein) and nu i ncrioal modcla (e .g., | K ft# 



isons of other cosm ological models). The observational data 
were collected from Gjjlego_et^jaLJ4^ IJUyet^ . (1996 



Connolly et al. (1997)t Bawicki et al. (1997)1; [Hughes et al 
1998)1; iGlazebrook et al. (1999)1 Madau et al. (1998)1; [Pas' 



19921; ^[ toinmota fo MiiUor 1005]| ^opoa ot oih 1007[). In thio 



section 
simulat 
a nd pIt 



irill moiko 



oatimato of E^T^Ei^.nmonoal ZSlt^^^l^^) tl'^l £.^.'!^i°'./l?5sF!!Tk "j^l 
Tho numorionl toohninuoa 



of gulp 



format 



(1998)1; ICowie et al. (1999)| ; ^teidel et al. (1999)| , and [Yan et 



sical ingredients of the model are described in Yepes 



et al. (1997) . The model includes a self-consistent treatment 
of the dark matter and baryonic components and effects of 
cooling, star formation, and SNe feedback. Simulations in- 
clude a multi-phase model of interstellar medium. Note that 
these simulations account only for SN type II so the con- 
tribution from SNI is therefore neglected in the estimate of 

The ideal simulation for our purpose would be a full 
modelling of cluster formation that would include forma- 
tion of the cluster galaxies, their starformation and feedback. 
However, with the numerical code used here, this would re- 
quire a significant sacrifice in the spatial dynamic range and 
mass resolution and would make it impossible to follow re- 
liably the starformation and feedback processes. Such simu- 
lation awaits future higher dynamic range simulations using 
adaptive mesh refinement technique. We choose the follow- 
ing compromise. We use many small-box galaxy formation 
simulations to determine statistically the number of super- 
novae, A^sN, that is expected to explode in a galaxy of a 
given absolute magnitude, Mb. We then use this relation 
and assume a galaxy luminosity function in clusters to esti- 
mate how many supernovae could have exploded in a cluster 
of a given mass. The energy released in these SN explosions 
would provide an upper limit on the amount of energy avail- 
able for IGM heating. While the number of SNe could be es- 
timated by assuming a particular A^sn — Mb relation, the use 
of simulations in this study spares us from making this ad- 
ditional assumption. As we describe below, the simulations 



al. (1999). All data points correspond to measurements of 
comoving UV or Ha luminosity densities. In order to trans- 
form to star formation densities, we have followed Madau's 



prescription (Madau et al. 1998) to correct the original data 
for dust extinction and to transform the luminosity densi- 
ties to star formation densities. All data points were prop- 
erly rescaled to the ACDM cosmological model used in the 
numerical simulations. The figure shows that the simula- 
tion results are in agreement (within the errors) with the 
observed evolution of the global star formation rate. The 
star formation rate in the simulations may actually be a lit- 
tle higher than the average observed rate, implying hence a 
larger number of exploded SNe. As we derived the SN rate 
from the galactic star formation rate shown in Fig. 1, this 
figure may serve as an illustration of how the SN explosion 
rate evolves with time. 

We have analyzed two additional simulations to assess 
the effects of resolution and box size. Particularly, the effects 
of resolution were checked by re-running one of the 5 Mpc 
simulations with 256'^ grid cells and particles (i.e., with eight 
times better mass resolution and twice the dynamic range). 
We have not found any significant changes in the global star 
formation rate or in the predicted number of supernovae 
(see below). To test the effects of the box size, we ran a 
simulation of 8.4/i"^ Mpc = 12 Mpc box using 300^ grid 
cells and particles, which gives the same resolution as the 
128'^ 5 Mpc runs but in a 2.4 times larger box. The results 
of this simulations are shown in Figs. |l| and ^ together with 
the results of other runs. The figures show that results of 
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■ Madau et al. 1998 
o Pascarelle et al. 1998 
Treyer et al. 1998 

• Cowie et al. 1999 

• Lilly et al. 1996 

□ Connolly et al. 1997 
a Sawicki et al. 1997 
A Steidel et al. 1999 

• Hughes et al. 1998 (SCUBA) 
^ Yan et al. 1999 (z=1.3) 

^ Glazebrook et al. 1998 (z = 0.88) 
^ Tresse &c Maddox 1998 (z = 0.2) 
^ Gallego et al. 1995 (z = 0) 



2 



3 



4 



5 



Figure 1. Comparison of star formation history in the Universe (points), estimated from the UV and Ha luminosity densities, corrected 
for dust extinction, with the star formation history obtained in the gasdynamical simulations used in this paper. The solid lines indicate 
the average and scatter of the star formation history in eleven 5 Mpc runs, while the thick solid line shows the star formation rate in 
12 Mpc run. All of the data points have been converted to the ACDM cosmology with Qg = 1 — = 0.35 and the Hubble constant of 
Ho = 70 km s-i Mpc-i. 



the large-box simulation are in agreement with results of 5 
Mpc runs. 

As we mentioned above, to estimate E^j^ expected from 
galaxies which end up in a cluster, we make use of the corre- 
lation between absolute magnitude of a galaxy at z=0 and 
number of type II SNe exploded in this galaxy through- 
out its evolution. The number of type II SNe exploded in a 
galaxy of absolute magnitude M, Nsn{M), is computed as 
the fraction of gas mass converted into stars of mass > IOMq 



divided by the IMF- weighted mean SNe mass. We use the 



Salpeter (1955) IMF, for which these numbers are 0.12 and 
22Mq (using lower and upper integration limits of O.IMq 
and 125Mq), respectively. 

Figure ^ shows z = correlation Nsn{Mb) for galaxies 
formed in the eleven 5 Mpc and one 12 Mpc ACDM runs. 
The correlation at z — can be well fitted by a linear fit 
(shown by solid line) of the form log(A'^jv) — a + hMs, 
with a — 1.411 and h — —0.344. Figure pi shows evolution 
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Figure 2. Correlation of tlie z = B-band absolute magnitude of simulated galaxies and number of SNe type II exploded in the galaxy 
throughout its history. The open circles represent simulated galaxies in 5 Mpc runs, triangles represent galaxies in the 12 Mpc run, and 
the solid line shows a linear fit to the correlation. Number of SNe, Nsn, is estimated as Ngi^ = 0.12Mt /22Mq, where M* is the mass of 
gas converted into stars, 0.12 is the fraction of the mass converted into stars with mass > IOMq, and 22M0 is the IMF- weighted mean 
supernova mass. Both numbers assume a Salpeter IMF. 



of this correlation with redshift. The figure shows that by 
z = 2 the number of exploded SNe is predicted to be « 
3 — 5 smaller than the number exploded by ^ = 0. Note that 
galaxies in our simulations are either isolated or are located 
in poor groups. It can be expected that formation of cluster 
galaxies occurs somewhat earlier than that of galaxies in 
poorer envi ronments (by about Az ~ 1, see, e.g., Gottlober 



et al. 2300) and the results for z > 2 should probably be 
interpreted as z > 3 — 3.5 instead. 



To estimate the supernovae energy input in a cluster of 
a g iven virial mass, Myjr, we convolve Nsn{Mb) fit with 
the ^chechter (1976) galaxy luminosity function 



(t){M) = 0.4 In 100. x 



l + a 



10' 



0.4{i\f,-Af). 



(3) 



where normalization parameter 0* is assumed to be equal 
to Avir times its field value. The parameter A^ir is the ex- 
pected virial overdensity in a given cosmological model and 
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Figure 3. Correlation of tlie z = B-band absolute magnitude of simulated galaxies and number of SNe type II exploded in the galaxy 
prior to a given epoch z (z = 0: upper left panel; z = 0.5 lower left panel; 2 = 1: upper right panel; 2 = 2; lower right panel. Symbol 
labeling is the same as in Fig. El 



is 334 for the ACDM model ado pted for our estimate (e.g., 
Lahav et al. 1991; Eke et al. 1996). The energy input is thus 



Elr^iM^.r) = E,A^,r ( Y^'-) / NsN{M)(P{M)dM; (4) 



where Eg = eiiEsNii is the energy input of a single super- 
nova e: 
of the c 



integration, and e/j, Esnii have the same meaning as in the 
previous section. 

The parameters of the luminosity function of galax- 
ies i n clusters appe ar to be similar to those of the field 
LF ( [Trentham l"998| ) and we will therefore neglect possible 
small differences between cluster and field LFs and cluster- 
to-cluster variations. We adopt parameters Mb = —19.5 
and a = —1.2 of the Schechter luminosity function consis- 
tent with recent measurements of B-band LF in the field (Da 



ix plosion, ckiM) is the field LF, R„ir is the virial radius Costa et al. 1994 ; Zucca et al. 1997) and in clusters (Tren 



luster, Mb and Mf are the bright and faint limits of tham 1998 ). The faint-end slope a is somewhat steeper than 
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in L Fs from most of other field surveys (e.g., Loveday et al 



1992 



N/ arzkeet al. 1994 



Lin et al. 1996). However, the steep 
value a = —1.2 better matches the LF of cluster galaxies 
and the faint end slope of the LF of the simulated galaxies 
1999b| ). Therefore, we adopt this value 



Elizondo et al. 



( 

in our analysis along with the normalization of the field LF 



I field 



0.02/i-='Mpc"^ (Zucca et al. 1997). This value may 



be u nce rtain by a factor of two (see Table 1 in Zucca et al 
1997 |). '^ he Egj^ estimate presented below is proportional to 
0, and can be simply rescaled for other values. We use the 



integration limits Mg — 



-22 and M| = 



-14. The results 



are insensitive to adopting a brighter Mj, or a fainter Mf. 
For consistency, we use en = 0.1 and Esnii ~ 1.2x lO^^ergs 
adopted in the previous section. 



4 RESULTS 

Figures ^ and ^ show results of the estimates described in 
the previous two sections. In Figure ^we compare estimated 
energy input from SNe with the thermal energy of the ICM 
gas. Top row of Fig. ^ shows Esn estimate using observed 
ICM metallicities and model A for SNII yields, while the 
bottom row shows the same estimate for yield model B (see 
Table |l|). 

The thermal energy of the gas is computed as 



Eth = 67r- 



linip 



Pg(r)T(r)r dr, 



(5) 



where k is the Boltzmann constant, rup is the mass of pro- 
ton, ^ — 0.6 is the assumed mean molecular weight of the 
ICM plasma, Pg{r) and T{r) are its radial density and tem- 
perature profiles. A density profile pg{r) is assumed, and 
the temperature profile is calculated from the equation of 
hydrostatic equilibrium. In Fig. ^ we assume that gas is 
dist ributed similarly to dark matter, and is described by 
the Navarro et al. (1997) (hereafter NFW), functional form, 

(6) 



Ps 



(r/r,)(l-f r-/r,)2 

with appropriate scaling of parameter c = Rvir/vs with clus- 
ter virial mass (NFW). The observed distribution of the ICM 
gas is more often described by the /3-profile: 

Po 



p{r 



(7) 



[l-H(r/rc)2]3f/2' 

where Vc is the core radius, and the parameter /3 controls the 
outer slope of the distribution. The thermal energy of gas 
distributed with the above density profile (for values of rc 
and f3 consistent with the observed range) is only ^ 10 — 20% 
higher than Eth of the NFW-distributed gas, the difference 



indistin guishable on the scale of Fig. \^ 



gradient is allowed. This implies that if metallicity gradients 
exist in clusters, simple galaxy formation models with the 
Salpeter IMF will have no difficulty in accounting for the 
observed amount of metals in clusters. 

Figure ^ also shows estimates for the cluster AWM7 for 
which a large-scale abundance gradient has been observed 
and its parameters measured (Ezawa et al. 1997). The es- 
timate was made using the observed metallicity gradient: 



Z{r) = Zo[l + {r/r 



-3/3 



with Zo = 0.59 and /3z = 0.8, 



where the core radius rc — 57.5h~^ kpc is equal to that of 
the gas distribution. The gas distributi on is described by 
the similar /3-profile with /3g = 0.58 (see Ezawa et al. 1997 



for details). We have made estimates for different gas frac- 
tions, fb, within the cluster virial radius but the results are 
only mildly sensitive to a particular value of /;,; the esti- 
mates shown in figs. ^ and |^ were made assuming ft — 0.2. 
The energy input estimate for AWM7 lies even lower than 
the solid lines at this mass because to calculate the latter 
we have assumed a core radius of rc — WOh~^ kpc, which 
results in a higher mass of metals and consequently larger 
estimated Esn- 

For reference, table |^ gives predictions for the num- 
bers of type II SNe in clusters of different masses based on 
our model (line 1), as well as the total masses of Fe and Si 
inferred from observations with assumptions of metallicity 
gradient and uniform metal distribution (lines 3-6). The ta- 
ble also gives the number of SNe required to deposit 1 keV 
per gas particle into the ICM (line 2), estimated assuming 
that average supernova deposits 1.2 x 10^" ergs. The mass 
of metals predicted in our model can be easily obtained by 
multiplying the number of SNe in line 1 of table by the 
corresponding mass- weighted yield given in Table ^. Note, 
however, that the predicted number refers to the type II SNe 
only, and the predicted mass of metals is thus only due to 
SNIIe. 

The question we would ultimately like to address is 
whether the energy input from SNe can noticeably affect 
the thermal state of the ICM. To answer this question, we 
need to know what energy input is needed to account for the 
observed properties of the ICM. It is not completely clear 
what energy is required. However, on theoretical grounds 
(Kaiser 1991; Evrard & Henry 1991) it is known that model 



predictions are in better agreement with the data when gas 
is assumed to be preheated (by some non-gravitational pro- 
cess) at an early moment. The preheating results in gas evo- 
lution corresponding to a higher adiabat, which affects the 
evolution of the accreted gas (in particular, some of the ac- 
creted gas may avoid being strongly shocked). 

Numerical simulations ( Metzler fc Evrard 1994|; 



Navarro et al. 1995; Mohr & Evrard 1997; Pen 1998) and 



semi-analytic models of clus t er evolution ( ^avaliere et al 



Fif urc \^ shows that expected energy input from SNe is. 



depending on the assumed uniformity of the ICM metallic- 
ity, ~ 5 — 10% of the gas thermal energy for poor clusters 
[M^ir ~ lO"/i"^M0) and <5% for rich clusters {M^ir ~ 
10^^ Mq). In case of the strong metallicity gradients 
these numbers are ~ 3% and ^ 1%, respectively. The es- 



19971; Fozzi fc Norman 1999t [Balogh et al. 19994 [Valageas 
Silk 199^ Wu et al. 1999D confirm that preheating result. 



in cluster properties that are more in accord with obser- 



vations. For example, to simulate SNe heating Pen (1998) 



preheats the gas in his gasdynamic simulations by injecting 
1 keV of energy per nucleon of gas, or, for plasma with pri- 
mordial comp osition, « 0.5keV per gas particle. Cavaliere 



timates from both Si and Fe agree very well between each 
other, for fpc ~ 0.5. Supernova energy input estimated from 
the numerical simulations agrees well with Esn derived from 
observed metallicities in the case where a strong metallicity 



et al. (1997) assume in their model that SNe preheat the 
intergalactic gas to temperatures of « 0.5 — 0.7 keV, which 
corresponds to (3/2)fc r ~ 0.75 — l.OkeV p er gas particle. 



Balogh et al. (1999a) and Wu et al. (1999) argue based on 
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Figure 4. Comparison of the energy input by SNe, i?5]v, and thermal energy of the ICM gas, -Ethi for clusters of different virial masses, 
Myir- The four panels show estimates of Egpf using observed metallicities of Si and Fe, as described in § ^; top row: using yield model 
A, bottom row: using yield model B (see Table Esn shown for the cases where the ICM metallicity is uniform throughout the 
cluster {dashed lines) and for the presence of a strong metallicity gradient {solid lines) (see § ^ for details). Solid triangles are estimates 
for the cluster AWM7, in which such gradient was observed. Esn was estimated for three different metal fractions contributed by type 
la SNe; = 0.0, 0.5, 1.0, corresponding to the shown groups of three curves and points (the top curves/points correspond to = 0.0 
and the bottom to = 1.0). Etfi{M^ir) {dotted line; the same in all panels) was calculated assuming the gas is distributed similarly 
to the dark matter (NFW distribution). E^i^ of the gas distributed with the "/3-model" is larger by ^ 10 — 20%. The energy released by 
type II SNe in the numerical simulations is shown by dot-dashed lines (the same in all panels). 



their semi-analytic calculations that the energy injection of 
~ 2 — 3 keV per particle is required to bring model predic- 
tions in accord with observations. 

We can compare our estimate of Esn to these numbers 
calculating the energy per gas particle as EsN/Np, where 



Np — fbM^ir / {i-imp) is the number of gas particles within 
the virial radius of cluster. Figure ^ shows EsN/Np for Esn 
estimated from Si and Fe and from galaxy formation sim- 
ulations. The figure shows that the maximum energy per 
particle of ~ O.lkeV can be injected by SNe if the ICM 
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Figure 5. Energy ner gas particle injected by SNe (= Egj^/Np, where A'p = fi,Myir/{fimp)). The curves and points have the same 
meaning as in Fig. U The upper row shows the estimate of Eg jv from the observed metaUicity assuming model A for SNII yields, while 
the bottom row shows the same estimate assuming model B. 



metaUicity is homogeneous, while in the case of a strong 
metaUicity gradient the typical energy per gas particle is 
only a few tens eV. In particular, the estimate of Esn for 
the cluster AWM7 is only ~ 0.002 - 0.01 keV per particle. 
The corresponding estimate from galaxy formation simula- 
tions is ~ 10~^keV per particle. These numbers are ~ 5 — 20 
times smaller than the typical energy injection assumed in 
the cluster formation models quoted above. 



5 DISCUSSION 

The results presented in § ^ allow us to assess the condi- 
tions required for the SNe energy input to be important 
in galaxy clusters. The primary conditions that are implied 
by our estimate of Esn from the observed abundances of 
Si and Fe are (i) large-scale uniformity of the metal abun- 
dances throughout the cluster volume and (ii) near 100% 
efficiency in transfer of the energy of SN explosion to the 
thermal energy of the IGM gas. The latter assumption is 
rather unlikely and the energies derived from the observed 
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Table 2. Estimates of numbers of SNe and masses in Fe and Si for clusters of different masses (all masses are in Mq). 



11 



Cluster mass 


1014 


5 


X 10^4 






1015 


Nsti (simulations) 


7.8 X 10» 


3.9 


X lO" 


7.8 


X 


lOlo 


^SN (foi' f ksV per particle) 


4.9 X IQii 


2.5 


X 10^2 


4.9 


X 


1012 


Mpe (observed, assuming Z-gradient) 


2.7 X 10^ 


8.i 


i X 10'' 


1.5 


X 


lOlO 


Mpe (observed, uniform distribution) 


1.4 X IQio 


7.0 


X lOio 


1.4 


X 


10" 


Mgi (observed, assuming Z-gradient) 


2.1 X 10^ 


6.9 X 10" 


1.2 


X 


lOlo 


Msi (observed, uniform distribution) 


1.1 X 10"' 


5.5 


X 10"' 


1.1 


X 


10" 



abundances should therefore be considered as the upper lim- 
its on the amount of SN energy that could have heated the 
IGM. 

There are but a few theoretical predictions and observa- 
tional data concerning the degree of uniformity of the metal 
distribution in clusters. Based on the numerical simula tions 
that in clude galaxy f eedb ack and metal enrichmen t, Met- 
zler & l^vrard (1994)1 and [Metzler fc Evrard (1997)| predict 



that large-scale metallicity gr adients should exist in clus- 

observed 



Ezawa et al. (1997) 



ters. On the observational side 
such a gradie nt in cluster AWM7. More recently, Finoguenov 
et al. (2000) reported similar large-scale (-R^O.5 — 1 Mpc) 
metallicity gradients detected using ASCA observations for 
several other clusters. It is not clear, however, whether such 



10 — 50 eV per gas particle, which is far short of the energy 
injection typically assumed to bring theoretical models in 
accord with observations: ~ 0.5 — 2 keV per particle. The 
estimate will still be short by a factor of ~ 5 — 10 even if 
100% of the energy of every SN explosion goes into heating 
the ICM gas. 

The above conclusions are for clusters of virial mass 
Af„ir> 10i4/i-i Mq. Figure | shows that the ratio of pre- 
dicted SNe energy input to the thermal energy of the ICM 
gas increases by about an order of magnitude as the mass 
is decreased from lO^^h'^ Mq to lO^/i'i Mq. This trend 
means that the SNe energy input may be much more im- 
portant for clusters of mass Mvir S 5 x 10^"^ Mq than for 

more massive cluster^. The mass 5 x IQi^/i^i Mq_ corre- 



gradien ts are ubiquitous. Our estimates of energy input from sponds to the ICM temperature of « 2 keV (e.g. 



observed metal abundances differ by a factor of 5 - 10 if al. 1996), while deviations from non-similarit y are observed 



we assume a uniform distribution of metals versus metal- 



Eke et 



licity g adients of the type observed in AWM7. New, deep 



observations of ICM metallicity profiles are therefore crucial 
to make this estimate much more reliable. With the launch 
of the Chandra X-ray satellite, such observations should be- 
come available. Our estimate of the SN energy input for 
AWM7 is two orders of magnitude lower than energy input 
which seem to be required to sufficiently preheat the ICM 
gas. 

Incidentally, the existence of large-scale abundance gra- 
dients in clusters would solve the prob lem of the total iron 

i pib 



mass in clusters. David (1997) and pibson et al. (1997) 



show 

that if the contribution of type I SNe to the iron production 
in clusters is relatively small, the total iron mass in the ICM 
is too large to be explained by type II SNe produced with 
Salpeter IMF. Brighenti & Mathews (1999) argue, however. 



that this solution is unattractive because it makes it difficult 
to explain the metallicities and radial abundance gradients 
in massive elliptical galaxies. It is clear from our analysis 
that the existence of large-scale abundance gradients in the 
ICM can reduce the estimate of the iron mass by up to an 
order of magnitude, thereby eliminating the need for a large 
number of SNII and flatter IMF. 

The predictions of the SNe energy input, Egj^, of the 
numerical simulations of galaxy formation presented in this 
paper, although consistent with observed evolution of the 
global starfomation rate in the Universe, are somewhat lower 
than the estimate from the metal abundances, Eg^pf. The es- 
timates Egj^ and iSJiv agree reasonably well if a metallicity 
gradient is assumed and the contribution of type la SNe to 
the iron enrichment is > 50%. In particular, the Egj^ esti- 
mate is actually higher than estimate of E^pf for AWM7. 
However, the energy input in this case is of the order of 



in re al clusters for tempera tures of < 2 keV ( Ponman et al 



199!;; Balogh et al. 1999a). Nevertheless, it appears that 



quantitatively our conclusions will stand for poor clusters. 
The entropy of the prehea ted gas required to e xplain obser- 
vations is ~ 100 keV cm ( [Ponman et al. 1999|| which corre- 
sponds to an energy of ~ 1.5(7ie/10~'^cm~'^)'^'''^ keV per par- 
ticle, where is electron number density. Thus, the energy 
injection into the gas in cluster cores (rie ~ 10"'' cm~^) is 
about 1.5 keV per pa rticle . Semi-analytical calculations of 
Balogh et al. (1999a)| and |Wu et al. (1999)| show that the 



energy injection required to explain the data may be even 
higher: ^2 — 3 keV per particle|^. 

Such energy input is marginally consistent with our 
E^Pf estimate in the case of uniform metallicities and e ~ 1. 
For the case of metallicity gradients, -EsV is more than an 
order of magnitude lower. The energy input, predicted 
from numerical simulations is even lower and is ~ 0.1 keV 
per particle even for e = 1. Therefore, the conclusion we 
draw from this analysis is that it is unlikely that the en- 



t Note that this conclusion depends on our assumption that to- 
tal luminosity of stars in clusters is proportional to the cluster 
mass (See § H). Although this assumption is reasonable, there is 
evidence that mass-to-light ratio of clusters and groups is a func- 
tion of system mass. The data indicates that mass-to-light ratio 
of galaxy groups is somewhat smaller than that of clusters (e.g., 
Bahcall et al. 1995). In this case our conclusion would not be 



changed. 



Balogh et al. (1999a) give an estimate of the required en- 

fcr/(Mmjfp2/3) 



tropy of preheated gas as S 



3.7 X 



10'^'' ergs g~5/3 cm^. This corresponds to Ep = LbfirriH p^^^ S ~ 
3.5 keV per particle if we assume a typical density of gas in cluster 
cores (p « 10~27 g cm~^) and the above value of entropy. 
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ergy in put from SNe is sufficient to preheat the intraclus- 



ter gas to the required entropy, unless all of the explosion 



ported by re cent direct numerical simulations of Thornton 
et al. (1998) who studied radiative losses of a SN remnant 



energy goes into heating of the gas and metal abundances 
are uniform throughout the ICM. Moreover, in light of the 
estimates, the SN energy input can only be impor- 



tant if starformation rate in cluster environments is a factor 
of 10 higher than th e average cosmic rate . Similar conclu- 
sions were r eached by Balogh et al. (1999a|) JValagcas&^Silk 
(1999)[]and |Wu et al. (1999)| . Recently, [Loowcnstcin (2000 



(SNR) for a grid of densities and metallicities of the ambient 
gas. The arguments and simulations, however, assume spher- 
ically symmetric evolution of SNRs in ambient gas of uni- 
form density. The efficiency may be higher if the topology of 
ambient gas density is very assymetric and the gas has been 
s wept up and preheated by previous, recently exploded SNe 
(Larson 1974). This, f or example, may be the case during a 



have also used observed abundance of Si in the ICM to esti- 
mate possible SNe heating and found that the implied SNe 
energies would not be sufficient to heat the entire cluster gas 
to the required levels (note that this estimate was done as- 
suming uniform distribution of Si and esN = 1). He pointed 
out, however, that SNe could still be the source of heating if 
only the gas in cluster cores was heated. In this case, heating 
would have to occur after or during formation of a cluster, 
not at early epochs as was assumed previously, but suffi- 
ciently early enough to be consistent with lac k of evolution 
of met al abundanc es at lower redshifts (2:^ 1; Mushotzky & 
Loewen stein 1997). Details and quantitative predictions of 
such a scenario are yet to be worked out. 



It i 3 obvious that there are a number of uncertainties in 



our estimates of the SNe energy input. The estimates of E^j^ 
made with the assumption of uniform ICM metallicity are 



strong starburst (e.g., Tenorio-Tagle fc Bodenheimer 1988 ). 
The parameter e is thus likely to be environment dependent 
and the average value would be determined by the relative 
number of SNe exploding during periods of quiescent star 
formation vs. the number of SNe exploding in starbursts. 
Regardless of the actual value, considerable radiation losses 
are expected and therefore it seems very unlikely that the 
efficiency is close to 100% (e = 1). 

Beside the problem of heating efficiency, it is also not 
clear how the heated interstellar gas and released SN en- 
ergy is transferred to the IGM (or ICM). Several transfer 
mechanisms have been suggested. Gas m ay be blown away 
frora galaxies by supernova- driven winds (Iviathews fc Baker 
197l|; lYahil fc Ostriker 1973t [Larson 1974| ) which subse- 



quently shock the IGM gas. Evidence for winds is indeed 
obser ved in some starburst galaxies (e.g., Heckman et al 



by a fac tor ~ 3 — 10 higher than the corresponding estimates 



in the case when a strong metallicity gradient is assumed. 



This uncertainty not only makes the -BJiv estimate uncer- 



199C ). However, only a small fraction of gas is expect ed to 
be blown away by sta, r bursts in massive galaxies (e.g., [Vlat 



teucci & Gibson 1995, Mac Low & Forrara 199!:) and there- 



tain, but also hinders comparisons of metal abundances pre- 
dicted by galaxy formation models with observations. This 
will likely be resolved in the near future with the advent of 
new, deep X-ray observations of clusters, but it is a major 
limitation at present. Currently, only one robust measure- 
ment of large-scale metallicity gradient has been obtained 
(Ezawa et al. 1997). This cluster, AWM7, conffims the ex- 



istence of strong metallicity gradients and the estimate of 
E^Pf for this particular cluster supports our conclusions. It 
is not clear, however, how universal such gradients are in 
clusters. 

Note also that our estimates are based on average abun- 
dances of Si and Fe from a large sample of clusters. Abun- 
dances in individual clusters may vary by a factor of 2 — 3. 
Thus, for example, abundances of Si and Fe (in solar units) 
vary in the range ■-^0.1 — 1 and ■-^ 0.15 — 0.45 , respectively 



( Mushotzky et al. 1996 Fukazawa et al. 1998). The energy 



fore ejected gas can only constitute a small fraction of ICM. 
Clearly, the same questions arise when we consider how en- 
ergy released by SNe can actually heat the IGM. If only a 
fraction of this energy is delivered to IGM, this effectively 
means a smaller value of esjv and strengthens our conclu- 
sions. 

The g as can also be transf erred to the ICM by ram pres- 
sure (e.g., Gunn fc Gott 197^ ) and tidal stripping. The effi- 
ciency of ram pressure in clusters is not well known. Recent 
numerical simulati ons, however, sug gest that it may actu- 
ally be rather low (Abadi et al. 199S;). The tidal stripping is 
probably the most efficient mechanism of delivering ISM gas 
to the intraclu ster medium, espec ially for low surface bright- 
ness galaxies (Moore et al. 1999). However, in this case the 
gas is transferred to the ICM relatively late, after the epoch 
of cluster formation, when a sufficiently deep potential well 
is formed. This is in conflict with high metal abundances ob- 



estimat es for individual clusters may therefore also vary by 



a corre! ponding factor. 



served in high-redshift cl usters (Mushotzky & Loewenstein 
1997|; iHattori et al. 1997|). 



The theoretical yields of Si and Fe from type la and 
type II SNe used in our analysis depend on specifics of the 
explosion model. T he Si yields from SNIa may be uncertain 
by a factor of two ( Nomoto et al. 19971 ; Nfagataki & Sate 
199 ^), j vhile all models predict similar (to ~ 10%) yields of 
iron. The yields of SNII for Si and Fe vary by ~ 30 - 40% 
between different models (e.g., Nagataki & Sato 1995). Yield 
models A and B used in our analysis approximately repre- 
sent the range of predictions and should therefore provide a 
fair estimate of uncertainty. Our conclusions hold for both 
yield models. 

The fraction of supernova explosion energy that can be 
available fo r gas heating is also rather uncertain. Larson 
(1974)[Jind Babul & Rees (1992) give analytical arguments 
that this fraction should be ~ 0.1. These arguments are sup- 



Recently, Snedin & Ostriker (1997) suggested that 
metal-enriched gas can be ejected at early epochs during 
galactic mergers. This mechanism may transfer metal-rich 
hot interstellar gas into IGM, where it can be further heated 
by shocks developed during a merger or after an encounter 
between ejected material and the ambient IGM gas. Despite 
the abundance of possible processes, it is not clear which 
process (or combination thereof) is responsible for the trans- 
fer of gas from galaxies into the intergalactic medium. It is 
clear, however, that this question needs to be clarified if SNe 
are to be considered a viable source of IGM heating. 

We have made a number of assumptions to estimate 
Egj^ from the galaxy formation simulations. Changing some 
of these assumptions can change the energy input estimate. 
First of all, our assumption of Salpeter IMF directly affects 
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the number of SNe per given mass of formed stars. IMFs 
Flatter than a Salpeter result in a larger number of super- 
novae and thus in a larger energy input for the same star 
formation rate. For instance, a 10% flatter slope with respect 
to Salpeter's results in a 50% increase in the number of SNe, 
given the same low-mass limit of the IMF. Indeed, a flatter 
IMF has been suggested as an ex planation fo r the observed 
iron abu ndances in cluster s (e.g. 



s tructed from massive ellipticals and small galaxy gr oups 
([David 1997|; [Renzini 199^; [Brighenti fc Mathews 1999|). El- 



David 1997 



Gibson et al 



1997|)~J[owever, note that [Brighenti fc Mathews (1999) 



ar- 



gue that flatter IMF is not consistent with the evolution of 
elliptical galaxies. The number of SNe depends also on the 
low-ma ss limit of the IMF, although in a less sensitive man- 



ner. Thus, an increase of the lower-mass limit by a factor of 
2 (from 0.1 to 0.2 Mq) results in an increase factor of 1.3 in 
the number of SNe. 

To calculate the number of SNe exploded during a Hub- 
ble time in all cluster galaxies we have assumed that the 
number density of galaxies in clusters is equal (1 -|- Avir) 
times its field value. This means that clusters represent the 
same fluctuation in number of galaxies as in their total mass. 
Although this is a reasonable assumption, we note that in 
the ACDM model (as well as in other low-matter density 
CDM cosmologies) studied here, a certain amount of anti- 
bias (6 ~ 0.5) is required for the model to be con s istent with 



observed galaxy clustering ( Kl ypin et al. 



1996 



Jenkins et 



al. 199^ Kravtsov fc Klypin 199£ ). This anti-bias arises pri 



r narily in the densest regio ns of galaxy groups and clusters 
(Kravtsov fc Klypin 1999). For an anti-bias of 6 ~ 0.5 the 



number density of galaxies would be two times lower than 
assumed in our analysis, which would reduce the estimated 
energy input by a factor of two. 

We neglected possible differences between the shape of 
the field and cluster luminosity functions. These differences 
a,ppear to be rath er small for the _B-magnitude LF used here 
( Trentham 1998 ), and we therefore think that the uncer- 
tainty associated with this assumption is relatively small. 

A more important assumption is that the global star for- 
mation histories of field and cluster galaxies are sim ilar. At 
present , there is no convincing evidence otherwise. Balogh 
et al. (: j999b)[ for example, argue that star formation activ- 
ity in cluster galaxies is not very different from that in the 
field. They argue, in fact, that field galaxies may produce 
more stars (and more type II SNe) than cluster galaxies in 
which the star formation is being gradually turned off after 
their infall onto cluster. This is in fact consistent with the- 



oretical predictions of Kobayashi et al. (1999) who present 



models for the evolution of the SN rate in clusters and the 
field. They predict that the rate in clusters is higher than in 
the field only at z^3.5, while at lower redshifts it is actu- 
ally lower due to a decreased contribution from SNe in spi- 
ral galaxies. Their predictions for the overall starformation 
rate in clusters are almost an order of magnitude lower than 
the starformation rate in the simulations presented here at 
z < 3.5 and are higher at higher redshifts. It seems unlikely, 
however, that their prediction can account for the required 
tenfold increase in number of SNe because only a small frac- 
tion of SNe in cluster galaxies explode at ^ 4. We therefore 
conclude that possible differences in starformation histories 
between cluster and field galaxies are too small to change 
our conclusions. 

Nevertheless, it is known that rich clusters have proper- 
ties different than if they would have simply had been con- 



lipticals and galaxy groups appear to have smaller gas frac- 
tions and lower metal abundances than rich clusters do. In 
particular, the ratio of iron mass in the ICM to the total 
blue luminosity of cluste r galaxies is c onsistently higher for 
clusters than for groups (Renzini 1997). It appears also that 



parameters of the models of elliptical galaxies that are tuned 
to produce the observed metal abundances in the ICM are 
inconsistent with abundance measurements in individual el- 
lipticals, the problem which cannot be solved by adjusting 



the SNIa contr ibution to the metal enrichment (Brighenti fc 
Mathews l"99s| ). These problems may indicate that an im- 



portant component is missing in our understanding of the 
ICM enrichment history and cluster evolution. However, our 
conclusions about the importance of SNe energy input can 
only change if the star formation rate in the volume from 
which the cluster forms is significantly higher at all epochs 
than that star formation rate in the field. 



6 CONCLUSIONS 

We have presented estimates of the possible energy input 
by supernovae into the intracluster medium. Although these 
estimates are prone to a number of uncertainties, we have 
defined conditions which determine whether SNe can be a 
significant source of ICM heating. The following main con- 
clusions can be drawn from our analysis. 

The SNe energy input, E^pf, estimated from observed 
ICM abundances of Si and Fe is only significant (~ 1 keV per 
particle) when we assumed that the distribution of metals 
in the ICM is uniform (no significant radial gradients) and 
that ~ 100% of individual SN explosion energy goes into 
heating the ambient gas followed by negligible cooling (e ~ 
1) (see § If large-scale metallicity gradients are assumed 
in clusters, the estimated energy input is ~ 0.1 — 0.5 keV 
per particle for e = 1 and, correspondingly, 0.01 — 0.08 keV 
per particle for a more realistic value of e — 0.1. 

As an example, we present estimates of the energy input 
for the cluster AWM7 for which the abundance gradient has 
been measured. We find that the observed abundance of iron 
in this cluster implies a SNe energy input of < 0.01 and < 0.1 
keV per particle for e — 0.1 and e — 1, respectively. 

The energy input, estimated using self-consistent 

three-dimensional numerical simulations of galaxy formation 
which include effects of shock heating, cooling, SN feedback, 
and multi-phase model of ISM, are ~ 0.01 and « 0.1 keV per 
gas particle for values of efficiency parameter e = 0.1 and 
e = 1, respectively. These values are somewhat lower than 
the values of E^pf (but are in good agreement with estimates 
for the AWM7). Nevertheless, the two estimates agree rea- 
sonably well if the existence of large-scale abundance gra- 
dients is assumed in clusters. We therefore emphasize the 
importance of new measurements of large-scale metallicity 
gradients for testing the theoretical models. 

Our estimates of the SN energy input in all cases, except 
the case of uniform ICM abundances and e = 1, fall short of 
the energy injection of ~ 0.5 — 3 keV per particle required to 
bring theoretical models of cluster formation in accord with 
observations. This suggests that supernovae are unlikely to 
be the only source of the ICM heating and should possi- 
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bly be supplemented (or substituted) by some other heating 
mechanism. Similar conclusions have been reached in recent 



stud ies of iBalogh et al. (1999al IValageas fc Silk (1999) 



and Wu et al. (1999), Valageas & Silk (1999) propose radia- 
tion from quasars as an alternative heating mechanism. This 
opens discussion of new possible processes for what appears 
to be a required high-redshift preheating of the intergalactic 
medium. 
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